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We investigate the influence of the appearance of excitonic states on van der Waals interactions 
among two Rydberg atoms. The atoms are assumed to be in different Rydberg states, e.g., in the 
|ns) and \np) states. The resonant dipole-dipole interactions yield symmetric and antisymmetric 
excitons, with energy splitting that give rise to new resonances as the atoms approach each other. 

Only far from these resonances the van der Waals coefficients, can be defined. We calculate 
the Ce coefficients for alkali atoms and present the results for lithium by applying perturbation 
theory. At short interatomic distances of several pm, we show that the widely used simple model 
of two-level systems for excitons in Rydberg atoms breaks down, and the correct representation 
implies multi-level atoms. Even though, at larger distances one can keep the two-level systems but 
in including van der Waals interactions among the atoms. 

PACS numbers: 32.80.Ee, 37.10.Jk, 71.35.-y 


I. INTRODUCTION 

Rydberg atoms, i.e. atoms with large principal quan¬ 
tum numbers n 3> 1, have been the subject of exten¬ 
sive research due to their unique optical and electrical 
features [l|. They are an ideal platform for the study 
of a wide range of quantum phenomena, mainly since 
their properties scale strongly with n, and often turn¬ 
ing them extraordinary relative to those of ground-state 
atoms. The property which has stimulated the most re¬ 
cent experimental and theoretical work is strong interac¬ 
tions among Rydberg atoms that are sraarated by rela¬ 
tively large distances of several pm [2-01 • 

On the other hand, electronic excitations can be de¬ 
localized among distant atoms through resonant dipole- 
dipole interactions to form collective electronic excita¬ 
tions that are termed excitons ii- Different processes 
of resonant energy transfer are possible m , here we con¬ 
centrate on the exchange type, e.g. ns + np O np + ns. 
Excitons have been introduced to a system of ultracold 
atoms in an optical lattice with lowest excited states 
where van der Waals interactions are negligible 0[ll. 
The formation of excitons in a cluster of Rydberg atoms 
has been investigated, but van der Waals forces were 
completely neglected and the discussion limited to two- 
level atoms |13l - [l^ . In previous work we investigated 
the influence of van der Waals interactions on the for¬ 
mation of excitons in an aggregate of two-level Rydberg 
atoms [l^ . Coherent energy transfer among Rydberg 
atoms that induces by resonant dipole-dipole interaction 
has been realized experimentally at large interatomic dis¬ 
tances of tenths pm where van der Waals forces are neg¬ 
ligible [13, [11. 

An excitonic state contains at least two different 
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atomic states. The simplest case is of two atoms in which 
one in the |ns) state and the other in the \np) state. 
The discussion can be limited to dipole-dipole interac¬ 
tions, where dipole-quadrupole, quadrupole-quadrupole 
and higher order interactions are neglected. Approximate 
long-range potentials can be derived by applying pertur¬ 
bation theory up to the second order in the dipole-dipole 
interactions. The lowest order term of the perturbation 
series results in resonant dipole-dipole interaction of the 
form where R is the interatomic distance with 

the resonant dipole-dipole coefficient . This interac¬ 
tion leads to a coherent mixing of the two possible states, 
which are \ns,np) and \np,ns). The second order term 
is of the van der Waals type of the form /R^, with 
the van der Waals coefficient Cg^. As was shown in our 
previous work [11, van der Waals interactions result in 
energy shifts that significantly influence the formation of 
excitons when atoms approach each other. Higher order 
terms play significant roles and can change completely 
the long-range interaction potentials [^ [^. Further¬ 
more, the appearance of resonances breaks down the va¬ 
lidity of the perturbative calculation and then other tech¬ 
niques are required, e.g. in using direct Hamiltonian di- 
agonalization [53]. Moreover, non-adiabatic interactions 
between Rydberg atoms in different electronic surfaces 
can be important pl . 

Long-range van der Waals interactions among pairs 
of Rydberg atoms have been intensively investigated, 
mainly using a perturbative approach [23l - l2a |. The dis¬ 
persion coefficients of the type Cl®, Cf^ and C|^, are cal¬ 
culated and listed for homonuclear dimers of alkali metal 
atoms in the ns — ns, np — np and nd — nd states, where 
both atoms are in the same state. These coefficients are 
of importance for experimental and theoretical applica¬ 
tions in strongly interacting Rydberg gases, especially 
for the current cold and ultracold Rydberg atom exper¬ 
iments, e.g., in the implementation of dipole blockade 
phenomena for quantum information processing 0,13, 
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But coefficients of the mixed type in which atoms are in 
different electronic states, e.g., Cg^, are not much em¬ 
phasized, while they are of importance for processes that 
involve resonant energy transfer [l^ . 

In the present paper we study the influence of the for¬ 
mation of excitons on van der Waals interactions among 
Rydberg atoms. We check the validity of using the simple 
model of two-level systems to describe excitons in inter¬ 
acting Rydberg atoms. We start by developing simple 
models using three and four atomic levels that provide 
a qualitative understanding of the different interactions 
among Rydberg atoms. In treating ns — ns and np — np 
Rydberg atoms we extract the limit in which atom-atom 
interactions have the forms /R^ and /R^. For the 

case of ns — np Rydberg atoms we derive the condition 
in which the interactions can be described by /R? 
and /R^ terms. Afterwards, we exploit perturbation 
theory in order to get quantitative values for all C 3 and 
Cq coefficients by summing over all the allowed atomic 
states. We treat Rydberg alkali metal atoms, and as an 
example we present the values for lithium atoms. 

We emphasize here van der Waals interactions of the 
mixed type with coefficient Cg^, and examine the effect 
of resonant dipole-dipole transfer on van der Waals in¬ 
teractions. The point is that as the two states \ns,np) 
and \np,ns) are degenerate, perturbation theory implies 
the removal of these degeneracy. But the states are cou¬ 
pled by resonant dipole-dipole interactions and the di- 
agonalization mixes and splits them to yield symmetric 
and antisymmetric orthogonal states. Then we use these 
orthogonal states as the zero order eigenstates in the per¬ 
turbative calculation of the van der Waals interactions. 
We found that as the symmetric-antisymmetric splitting 
energy is i?-dependent new resonances appear as the in¬ 
teratomic distance decreases. 

The paper is organized as follows. In section 2 we 
present a qualitative study of van der Waals interactions 
among two Rydberg atoms using multi-level simple mod¬ 
els. Quantitative derivations of van der Waals interac¬ 
tions appear in section 3 using perturbation theory. Sec¬ 
tion 4 contains calculations of van der Waals coefficients 
for alkali Rydberg atoms, and the results are presented 
for lithium atoms in section 5. A summary is given in 
section 6 . The appendix includes the angular parts of the 
dipole moment matrix elements. 


II. RYDBERG ATOM INTERACTIONS: SIMPLE 
MODELS 

We start by treating Rydberg atoms using models of 
two, three and four-level systems. The derivations pro¬ 
vide a qualitative understanding of the type of interac¬ 
tions that can appear among Rydberg atoms and pave 
the way towards the quantitative treatment presented in 
the next section. We examine the validity of using two- 
level systems for describing excitons in Rydberg atoms. 
The model of two-level systems is widely used to describe 


Frenkel excitons in organic solids involving lowest ex¬ 
cited states where van der Waals interactions result in 
small energy shifts ii- But for Rydberg atoms van 
der Waals forces are significant and the appearance of 
resonances is probable. Hence we examine the validity of 
using the simple model of two-level systems for Rydberg 
atoms. Our objective in this section is to extend the sim¬ 
ple model into multi-level systems, which are necessary in 
order to exploit the appearance of resonances. We show 
that far from resonances one can keep the simple model 
but in including energy shifts due to van der Waals in¬ 
teractions, which implies quantitative calculations of the 
van der Waals coefficients that we calculate in the next 
section using perturbation theory. 

A. Two atoms in the same Rydberg state 

Let us assume the two atoms to be in the same internal 
atomic level, say in the (nl) state. The two-atom state is 
Itt) = \nl,nl) of energy Ei = 2Eni. We assume a single 
channel to be close to resonance with this state. Namely, 
we consider a process of the transfer type nl + nl o 
n'l' + n”l". We have two possible degenerate final states 
\pi) = \n'V,n"l") and \p 2 ) = n'/'), with the energy 

Ef = En'i' +En"i"- The states have the energy detuning 
A = Ef — Ei = En'i' + En"i" — ‘2Eni. We assume here 
coupling among the [tt) state and each one of the \pi) 
and IP 2 ) states, with the coupling parameter J, which we 
specify later. We neglect coupling among the \pi) and 
IP 2 ) states. The Hamiltonian that is restricted to the 
states defined above can be written as 

H = Ei\-K){'K\+ Ef {\pi){pi\ + \p2){p2\) 

+ + k)(P2| + |pi)(7r| -h |p2)(7r|) . (1) 

In matrix elements we have 

/O J J\ 

H = Ed+\JA0\. (2) 

\J 0 A J 

We diagonalize the matrix to get the characteristic equa¬ 
tion (A — A) (A^ — AA — 2J^) = 0, with the three solu¬ 
tions 

Ai = ^-ikA2 + 8J2, A 2 = ^+ix/A2 + 8J2, A 3 = A. 

(3) 

In the far off-resonant limit, that is A ^ J, we get 

Ei!=sE,-2^, E 2 -Ef + 2^, E 3 = Ef. (4) 

Transfer processes are induced by resonant dipole-dipole 
interaction that have the form J = — , where R is the 

interatomic distance, and C 3 will be calculated in details 
later. We can also define D = .ffi, which is of 

the van der Waals type. In the general case we have D = 
— where one need to calculate the van der Waals 
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FIG. 1: The three states |7r), \pi), and |p 2 ) are illustrated, 
with the coupling among them. 


coefficient Cq for each specific states by considering all 
possible atomic transitions, which is the main aim in the 
next sections. Next we consider the case of ins) state in 
details. 


In the case of [tt) = \ns,ns) of energy Ei = 2Ens, we 
have, e.g., the channel ns ns O np + (n — l)p, that 
is |pi) = \np,{n — l)p) and \p 2 ) = \{n — l)p,np), with 
the energy Ef = Enp + and with the detuning 

Ass = Enp + E(^n-i)p - (see Fig. [T]). The energy 

of the two interacting atoms is Ess = 2Ens + Dss , with 
the van der Waals interaction Dss = —2 ^^^ ^ In 


the general case we can write Dss = Later on we 

calculate the Cl® and C|® coefficients in details. Similar 
consideration holds for the case of two (np)-state atoms. 



FIG. 2: The four states |7ri), |7r2), |pi), and |p 2 ) are illus¬ 
trated, with the coupling among them. 


In matrix elements we get 


H = Eii + 


/ 0 J J' 0 \ 

J 0 0 J' I 

J' 0 Asp 0 

V 0 J' 0 Asp / 


( 6 ) 


We diagonalize the Hamiltonian, to get the characteristic 
equations 


B. Two atoms in different Rydberg states 


(A - Asp)" (A^ - J2) _ 2A (A - Asp) J' 2 + J' 4 ^ 0, (7) 


Here we consider the case in which the atoms are in 
different states, e.g., one atom in the (ns)-state and the 
other in (np)-state. We get two degenerate states Itti) = 
\np, ns) and |7r2) = |ns, np) of energy Ei = Ens+Enp. We 
have, e.g., the possible transfer channel of ns+np o (n — 
l)p -I- (n -I- l)s. Further, we consider the two degenerate 
states \pi) = |(n-|-l)s, (n—l)p) and |p 2 ) = \{n — l)p, {n + 
l)s), of energy Ef = E(^n+i)s + L^(n-i)p, (see Fig. [2|). 
Their energy detuning compared to [tt) state is Asp = 
E(^n+i)s-Ens + E(^n-i)p-Enp- The states Itti) = \np,ns) 
and |7r2) = \ns,np) are coupled by the resonant dipole- 

dipole interaction of strength J, where J = —The 
state Itti) = \np,ns) is coupled to \pi) = \{n + l)s, (n — 
l)p) by the resonant dipole-dipole interaction of strength 
J', and the state |7r2) = \ns,np) is coupled to \p 2 ) = 
|(n — l)p, (n -I- l)s) with the same parameter J'. We 
neglect the resonant dipole-dipole interaction among the 
|Pi) = 1(41 + l)s, in-l)p) and |p 2 ) = |(n- l)p, {n+l)s) 
states, due to their small amplitudes. 

The Hamiltonian can be then written as 


with the four solutions 



A ^SP ^ J 1 

=2 2^ 

/(Asp± J)"+4J'2, 


A3,4 = 2 +2^ 

/(Asp± J)"+4J'2. 

(8) 


In the far off-resonant limit, that is Asp ^ J', we get the 
energies 

J' 2 

Ea.b — Ejis + Enp ^ J “T j-i 

^sp ^ •J 

J' 2 

Ec4 = E(^n+i)s + C(n-i)p + “T—v-r- (9) 

^sp ^ ^ 

As far as Asp ^ J the shifts are of the van der Waals 
type, where Dgp = — J' '^/Asp. In general we can write 

Dsp = —flf -1 where the Cq coefficient will be calculated 
later. 


C. Two-level atoms 


H = Ei[\TTi){-Kx\ + \'n2){-K2\) + Ef {\pi){pi\ + \p2){p2\) 

+ J (ki)(7r2| + |7r2)(7ri|) 

+ (ki)(pi| + k2)(p2| + ki)(7ri| + |p2)(7r2|) . (5) 


In the light of this result, we can go one step back 
and start with two effective two-level atoms including 
van der Waals interactions. We again treat the case of 
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FIG. 3: The two states |7ri) and \'ii 2 ) are illustrated with the 
coupling J among them. 


reads H = Hi-\- H 2 + V, where Ha is the a independent 
atom Hamiltonian, and V is the interaction among the 
two atoms. We assume the independent atom eigenstates 
to be known and given by where the 

eigenstates are orthonormal with {4>^^\4>a) = SapSij, and 

The lowest order interaction among neutral atoms is 
the dipole-dipole one, where 


one atom in the (ns) Rydberg state, and the other in the 
(np) state, as seen in Fig. |31 The atoms are separated 
by the distance R. We have two states [tti) = |np, ns) 
and \'K 2 ) = \ns^np), and they are degenerate with the 
energy Eq = Ens + Enp + Dgp, where we include the 
van der Waals interaction among the two atoms D^p = 
—hC^/R^. The energy transfer parameter among the 
two states is as before J = —hC^^/R^. These consider¬ 
ation exactly fits with the above four-level model in the 
limit of Asp ^ J, where the detuning Agp is much larger 
than the symmetric-antisymmetric splitting. 

The Hamiltonian is now written as 


V = 


1 


47reoi?^ 


(Ai ■ A 2 ) 3 


(Ai ■ R) (A2 ■ R) 

i?5 


(13) 


with R = Ri — R 2 , and the interatomic distance is i? = 
|R|. Here fia is the dipole operator at atom a. We 
assume no permanent dipoles to exist, and hence /x” = 
(Klfj-JK) = 0- We have with (i ^ j). 

The atoms are taken to be localized along the 2 -axis, 
which is also the quantization axis, hence 

V = + AiyA2y - ■ (14) 


H = Eo (|7ri)(7ri| -|- |7r2)(7r2|) -b J (|7ri)(7r2| -b |7r2)(7ri|). 

( 10 ) 

The Hamiltonian can be diagonalized by using the col¬ 
lective states 

I , > kl) -b |7r2) ki) - k2) 

ka) = ——, kb) = —, (11) 

which gives 

H = Ea ka)('0a| + \tjjb) {tpb\, (12) 

with the energies E^ = Eq + J and Eb = Eq — J, which 
fit exactly with the above derived energies. 

As other states can be close to resonance with the con¬ 
sidered states, the picture of two-level atoms breaks down 
at short interatomic distances of several fj,m. Then, the 
excitonic picture implies multi-level systems, as beside 
the considered states, e.g. ns and np, the formalism 
must include all close to resonance states, e.g (n -b l)s 
and {n — l)p, as treated before. The model of two-level 
systems can be reserved for larger interatomic distances 
when other states are off-resonance, but in including van 
der Waals interactions among the atoms. We adopted 
this direction in treating N Rydberg atoms in our previ¬ 
ous paper [l^. At much larger interatomic distances of 
tenth of fim the van der Waals interactions are negligible. 

After this qualitative study of the resonant dipole- 
dipole and van der Waals interactions, we move to quanti¬ 
tative calculations of the different C 3 and Cq coefficients. 


We apply perturbation theory in order to calculate the 
resonant dipole-dipole and van der Waals coefficients by 
deriving the corrections to each dimer energy. 

We write the perturbative terms up to the second or¬ 
der explicitly where the two atoms are assumed to be in 
the states i and j. The lowest order correction to the 
free energy is = (c))®, ^^ 2)1 where for the 

dimer we use a product basis ki,? 5 < 2 ) = kDk^)- This 
term give rise to resonant dipole-dipole interaction that 
is responsible for the energy transfer among atoms at 
different internal states. When the two atoms are in the 
same i-th quantum state one has {(j)\,(j) 2 \V\(j)\,(j) 2 ) = 0 . 
Only terms of (i k k) and (j k 0 nonzero, with the 

appropriate selection rules. We define ^—, 

where now 

~ , (15) 

and the transition dipole matrix element is defined by 
= ('('akarlki), with {r=X,y,z). 

Resonances can appear among the dimer states, e.g., 
for kk '(' 2 ) kik 2 ) states with the energy E 1 +E 2 = 
El + E 2 , which implies the use of degenerate perturba¬ 
tion theory. In the following we present separately the 
degenerate and non-degenerate perturbation theory. 

A. Non-degenerate states 


III. DIMERIC ENERGIES IN PERTURBATION 
THEORY 

We give first a general presentation of a perturbative 
treatment for the electrostatic interactions among Ryd¬ 
berg atoms. For a system of two atoms the Hamiltonian 


For non degenerate states kk '(' 2 ) l^e first order term 
(kk '/' 2 kkik 2 ) vanishes, because the atoms do not have 
a permanent dipole. The second order term is 

=-A- E^-Ei + E‘-Ei -' <“> 




















5 


which results in van der Waals interactions. Since V has 
a 1 /i?^ distance dependence, we can write 
where we have defined 


function. In general the dipole moment matrix element 
can be factorized as where 

the radial integral is 




( ’ IV 

\lEtx9^ + 

2 

II 

H 

Weo) ^ 

E^ -Ei+E!2-E^2 

(17) 

^0 

and the angular part is 


dr T R^i’i’Rjii^ 


( 21 ) 


The present perturbation theory can break down due to 
the appearance of resonances at + = E 2 + E\, 

and one appeals to other methods, e.g., to the direct 
Hamiltonian diagonalization M- 


-Im.l'm' 



d(j)sine Yi*^{e,(j))eiYi>jn'{d,(j)). 

( 22 ) 


Here are the axis unit vectors with 


B. Degenerate States 


We now consider states which a state 

with the same energy exists, i.e. Ei + Ej = 
Ek+Ei, and for which also = ((^^, ^ 0 . 

Then we have to use degenerate perturbation theory and 
first diagonalize the degenerate subspace. Here we con¬ 
centrate in degenerate states of the type 14 ) 1 , 4 ^) and 
|(/)j,(/)|). The ’correct’ zero order states have the form 

\4^t) = \4^\.4^i)±W,<4i), ( 18 ) 

with energy 

E^^=El+Ei±rE (19) 

The first order contributions then vanish and the second 
order term, which has a 1 /i?® behavior is calculated by 


D 


ij _ 
± — 


tr E^ + Ei-Et ’ 


( 20 ) 


Ga; = sin0cos((), Gy = sin0sin((), ez=cos9. (23) 

The radial wave function for alkali atoms in atomic 
units is given by 


Rni{r) = 


n 


2 I {n* - I* - 1 )! 


(n* + /*)! 


X 



e-r/u* 


r2/* + l 



(24) 


where n* = n — Sni and I* = I — 6ni + d(/), with the 
quantum defect Sni = ^0 + and I{1) is the near¬ 

est integer below or equal to Sni ■ For the energy we have 
Eni = 2 {n-s i)‘^ ■ Moreover we neglect fine and hyper- 
fine splitting, which is a good approximation for lithium 
atoms that we present later. 


B. C'l’’ and C|^ Coefficients 


where the sum over k and I includes all states but 
Note that resonances appear at E^ + E 2 = E'^ and then 
a special treatment is required. 


IV. INTERACTIONS BETWEEN ALKALI 
RYDBERG ATOMS 

We calculate the Cz and Cq coefficients for specific 
cases that include the and \4Y'^) states, and we 

concentrate in interactions among alkali Rydberg atoms. 
We use the nondegenerate perturbation theory in order to 
calculate C'|'’ and Cg^ coefficients. Next we calculate Cg^ 
coefficients, and as here a resonance appears between the 
states \4)Y i4h^) and |((>"^, (4^) we use degenerate pertur¬ 
bation theory. Let us first present the atomic states for 
alkali Rydberg atoms. 


A. Atomic eigenstates 

The atomic state is given by the wave function 
'fpnim{r,d,(4) = Rni{r)Yimld,(4), where Rniijr) is the ra¬ 
dial wavefunction and Yim{9, 4’) is the spherical harmonic 


We present here detailed calculations of C|® and 
coefficients for atoms in the (ns) and {np) states. We 
concentrate on alkali metal atoms excited to Rydberg 
states and hence we deal with a single electronic state at 
each atom that is represented by three quantum numbers 
(n, 1 , m), where n = 1 , 2 , • • •, Z = 0 , • • • , n — 1, and m = 

—I, • • • , +1. 

We write the expression for the van de Waals coeffi¬ 
cients as follows 




^np7n,npm 

where we defined 




Aupm 
■^npfh 1 


(25) 


A' 


nl_m 

nlfh 


= E 

n'V,n"l" 


Im^l' 


I'Tnl 'T' 


nl 


A 


nl,nl 
n'l' ,n"l" 


2 


(26) 


and the prime over the summation indicates that nl ^ 
n'l' and nl 7 ^ n"l". We have 


A nl,nl 


^n'V H“ Ejiffiff Ejii 


(27) 
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and 


Ip 


lm,l' 


\t^x 

m' ,m" 

, ,, ,, 2 
m ~lm,l m 

f^z 




+ K 


(28) 


As we treat two identical atoms, we dropped the atom 
index. We calculate the radial parts, I, of the A pa¬ 
rameters numerically, and the angular parts, |p|^, are 
calculated analytically and presented in the appendix. 

We aim now to calculate the Cg coefficient, which can 
be written as 

C'r™ = ^P^’"™Kr^ (29) 

where 


^s,pm ^ ^s,pm~pm,s ~s,pm^pm,s _ , (30) 

and here we have = — | and = i. 


C. Cg^ Coefficients 



FIG. 4: C|® as a function of n for lithium atoms. 


The i? dependence in the denominator through 
J®’^™(i?) can bring new resonances as atoms approach 
each other and then the present perturbation theory 
breaks down. Event hough, far from resonance, that is 
^ J®’P'"(i?), the coefficients can be safely 
defined by Z?®^ = —hC^/R^, where 


One needs to be careful when calculating the Cg^ co¬ 
efficients, as the states |ns,np) and \np,ns) are degener¬ 
ate, and perturbation theory implies the removal of this 
degeneracy, as we presented before. The two states are 
coupled by the resonant dipole-dipole interaction, which 
yields two diagonal states of symmetric and antisymmet¬ 
ric mixing, where these orthogonal states used as the zero 
order basis for the second order perturbation theory. 

The diagonal eigenstates are 


. , \ns,npni) ± |npm,ns) 

' ^ " V2 

(31) 

with the diagonal energies 


E^{R) = Er^s + Er^p ± J"®’"P’"(i?), 

(32) 

where 


1 ^ns,npm 

j"®’"p'"(i?) = % ■ 

(33) 

The diagonal states have splitting energy of 2 J"® 

,npm 


We aim to calculate 


T~\ns ,npm 
^Oi 


1 


47reoZ?^ 


X E 

k,l 


\.,ak,,(y.l _i_ , ,ctk , ,ocl _ fy,.ctk,,ocl\ 

iH'x f^x H'z 


Ek + El- E^{R) 


-,(34) 


with (a = ±), {k = n'l’m') and [I = n''l"m''). 
The calculation yields 


( Y 

V47reoi?3y 


E 

n'l' ,n"l" 


sX 

^pTn,i 


I'Tns 7 "' 


np I 
n"l"\ 


A ns,np 




(35) 


As before the radial parts are calculated numerically, and 
the angular parts calculated analytically and presented 
in the appendix. 


V. VALUES FOR LITHIUM ATOM 


We calculate now all Cg and Cg coefficients for lithium 
atoms, by using the known defect parameters [27|. For 
ns we have 6 g = 0.399 and 62 = 0.03, for np we have 
Sq = 0.0473 and = —0.026, and for nd we have Sg = 
0.002 and S!^ = -0.015. 


Converting the above A parameters to SI units requires 
the change A —>• A, as we have a factor of 

from the energy, and a factor of ag from the radial func¬ 
tion I. We obtain Cg""’”" = BA"® and = 

where B = Using eg ^ 0.0055 e/(yA), 

and oq ~ 0.53 A, we get B fa 0.6 eUA®, or B = hB, with 
h fa 4.135X 10“^® eVS, we get B fa 1.5x 10“^° Hz{pm)^. 
Moreover, we have where F = 

1^. We use F = hF with F fa 10^ Hz{fini)^ 

We plot now the results for all C|®, Cg^ and Cg^ co¬ 
efficients ranging from n = 50 up to n = 100. In Fig. 0] 
we plot C|®, and in Fig. [5] we plot Cg°’^°, and 

^p±i.p±i jjgj.g (jss coefficients are negative that lead to 
repulsive van der Waals forces, and Cg^ are positive that 
lead to attractive ones. The results for C|’^° and Cg^^^ 
are plotted in Fig. [51 Cg^^ coefficients are positive that 
lead to attractive resonant dipole-dipole interactions, and 
are negative that lead to repulsive ones. 
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50 60 70 80 90 100 

n 

FIG. 5: and as a function of n for 

lithium atoms. 



FIG. 6: and Cg’^^^ as a function of n for lithium atoms. 

Now we present the results for energies that de¬ 
rived in using the diagonal states as zero order in per¬ 
turbation theory for lithium atoms. The main feature of 
is in their R dependent. We present the plots for 
different interatomic distances. In Fig. [7] we plot 
and for R = 7fim, and in Fig. |8]for i? = 2/im. It 

is clear that resonances start to appear at smaller n as 
atoms approach each other. The and energies 
are different due to the symmetric-antisymmetric energy 
splitting, but for small n as the splitting is much smaller 
than the energy detuning then the coefficients become 
equal, and in this limit one can define Cg^ coefficients. 
In Fig. [9] we plot and coefficients from 

(n = 20) up to (n = 70). In this region no resonances 
appear down to small interatomic distances of 2 fim, as 
it is clear from the above plots. Here the C'g’^*'°^ coeffi¬ 
cients are negative for n > 30 and then the van der Waals 
interaction is repulsive, but the coefficients are 

positive and the interaction is attractive. 


VI. CONCLUSIONS 

In summary we investigated van der Waals interac¬ 
tions among two Rydberg atoms for the case where the 
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FIG. 7: and as a function of n for R = 7^tm of 

lithium atoms. 



FIG. 8: and as a function of n for R = 2^m of 

lithium atoms. 


atoms are in different internal states. For example, we 
treated in details the case with one atom in the (jis)- 
state and the other in (np)-state. The simplest model 
for achieving a qualitative understanding of the interac¬ 
tions is of four-level atoms, e.g., with |ns), \np), \n'p) 
and |n"s) states of energies Ens, Enp, En’p and En"s, 
respectively. The two states \ns,np) and \np,ns) are 
degenerate and couple by resonant dipole-dipole inter¬ 
action of the type Jsp = —hC^^/R^, which mixes and 
splits the states to give symmetric and antisymmet¬ 
ric orthogonal ones. We show that van der Waals in¬ 
teraction among the two atoms can be formulated as 
Dsp = -hCl^/R^ in the limit of off-resonance, that is 
in the limit \Ens + Enp — En'p — F^n"s| ^ I'^spl- In this 
limit the model of two-level atoms is valid, but in includ¬ 
ing van der Waals interactions. 

Next we presented a quantitative calculation of van 
der Waals coefficients among two alkali Rydberg atoms 
in considering all allowed electronic states. Due to 
the resonance among the states \ns,np) and |np, ns) 
we applied degenerate perturbation theory. The zero 
order states used in perturbation theory are the or¬ 
thogonal states \nsp)± = with energies 

j^nsp _ found that due to the 

symmetric-antisymmetric energy splitting dependence on 







































FIG. 9: and as a function of n for lithium 

atoms. 

R new resonances appear when atoms approach each 
other, and only far from these resonance one can define 
Cg^ coefficients. Moreover, using non-degenerate pertur¬ 
bation theory we calculated all van der Waals coefficients 
among alkali Rydberg atoms in the same state, mainly 
for the |ns, ns) and |np, np) states, and we presented the 
results for lithium atoms. 
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Appendix A: Angular Part of Dipole Moment 
Matrix Elements 


In this appendix we present the angular part calcu¬ 
lations of the dipole moment matrix elements. Direct 
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